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Abstract 

A class of nonlocal nonlinear wave equation arises from the modeling of a 
one dimensional motion in a nonlinearly, nonlocally elastic medium. The 
equation involves a kernel function with nonnegative Fourier transform. We 
discretize the equation by using Fourier spectral method in space and we 
prove the convergence of the semidiscrete scheme. We then use a fully- 
discrete scheme, that couples Fourier pseudo-spectral method in space and 
4th order Runge-Kutta in time, to observe the effect of the kernel function 
on solutions. To generate solitary wave solutions numerically, we use the 
Petviashvili’s iteration method. 
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1. Introduction 

In this article, we study the nonlocal nonlinear wave equation 

utt = (3*{u + g{u))^^, (1.1) 

where is a nonlinear function of u. Here * denotes the convolution in spatial 
domain and the kernel function (3{x) is an integrable function. Equation fll.ip 
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has been derived to model one-dimensional motion in an infinite medium 
with nonlinear and nonlocal elastic properties from nonlocal elasticity theory 
[^, 1^. Duruk et. al. have proved the global existence of the Cauchy problem 
corresponding to fll.ip together with initial conditions 

■u(a;, 0) = (/)(x), Ut{x,Q) ='ip{x) (1.2) 

with enough smoothness on the initial data, and continuity condition on the 
nonlinear function in j^. They assume that the Fourier transform of the 
kernel function (3{x) satisfies 

0 < ^{k) < ^(l + for all A; e R (1.3) 

for some constant C > 0, r G R and r > 2. Suppose that u and all its 
derivatives converge to zero sufficiently and rapidly as x —)■ ±cxo. For the so¬ 
lutions of the eq. flLlI) subjected to these boundary conditions the conserved 
quantity (energy) |3[ is given by 

E (t) = \\Put\f + \\u\f + 2 f G{u)dx. 

Jr 

Here the operator P is defined as 

\\k\Jmj 

r 

where P ^ is the inverse Fourier transform and G{u) = / g{s)ds. 

Jo 

For some suitable choice of the kernel function /5(x) the equation (II.Ih be¬ 
comes well-known nonlinear wave equations. In the case of the Dirac delta 
function, (3 = 6, the condition fll.3p is satisfied by r = 0 and the equation 
(mil is the well-known nonlinear wave equation 


Pt ^xx 9 ■ 

Choosing the exponential kernel j^, (3{x) = the equation fll.ip be¬ 

comes the improved Boussinesq equation (IBq), 

^XX ^xxtt {.divJ^XX • 
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Since j3{k) = (1 + k'^) we have r = 2. A final example is the double¬ 
exponential kernel 

Z[Ci C2) 

where Ci and C 2 are real and positive constants. As (3{k) = {l + rjik'^+ ri 2 k'^)~^ 
where tji = c\ + and r ]2 = cfc|, we have r = 4. In this case the equation 
fll.ll) becomes the higher-order Boussinesq equation (HBq), 

'^xx Vl^xxtt A ^2^xxxxtt (5^(^))a;a; • 


For these special cases there are numerous studies in the literature both 
numerical and analytical (see [ij and 


16 and the references there in). The 


existence and stability of solitary wave solutions of the nonlocal nonlinear 
wave equation fll.ip has been shown in a recent study j^. The nonlocal 
term makes it harder to investigate the equation analytically. Therefore the 
numerical solution is important to understand the evolution of solutions in 
time and the effect of the kernel on the solutions. 

The existence of the convolution integral and the condition fll.Sp on the 
Fourier coefficients of the kernel function /9(x) make the spectral methods a 
natural choice for the numerical solution. Spectral methods are widely used 
for local equations. Some of the numerical studies using spectral methods for 
the nonlocal equations are as follows: In BQ , the authors have studied a 
class of nonlinear nonlocal dispersive wave equation by using the combination 
of Fourier Galerkin spectral method and explicit leap-frog scheme. In [sl, the 
author has studied regularized Benjamin-Ono equation by using both Fourier 
Galerkin and Fourier collocation methods. 

In the recent studies UM a Fourier pseudo-spectral method has been 
used for the numerical solutions of the IBq and HBq equations. Thus the 
efficiency of the method has been tested for the special cases of this class of 
nonlocal equations. Therefore in this study we use Fourier pseudo-spectral 
method, together with a fourth order Runge-Kutta method, to solve fll.ip nu¬ 
merically. In Section 2, we give notations and preliminaries. In Section 3, we 
prove the convergence of the semi-discrete scheme. In Section 4, we introduce 
a fully discrete scheme combining the Fourier pseudo-spectral method and 
fourth order Runge-Kutta method for space and time discretization, respec¬ 
tively. Although the exact forms of solitary wave solutions are known for the 
IBq and HBq equations, we do not have the exact form of the solitary wave 
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solutions for the general nonlocal nonlinear wave equation. Therefore, we 
use the Petviashvili’s iteration method to construct the solitary wave prohle 
numerically in Section 5. Finally, we present some numerical implementation 
of the method in Section 6. 


2. Notations and preliminaries 

We use the function space and the Sobolev space H^. (•, •) denotes 

the inner product and || ■ || is the corresponding norm dehned, respectively 
as 

{u,v) = / u{x)v{x)dx, and ||n||^ = (n,n) (2.4) 

Jq 

where hi = (—T, L). The periodic Sobolev space is equipped with the 

norm 


Ml = 


^ (1 + k^)‘\utf 


(2.5) 


k=—oo 


The Banach space Xg = C^([0, T]; iJp(r2)) is the space of all continuous 
functions in Hp{Q) whose distributional derivative is also in Hp{Q), with 
norm ||n||x^ = rnax (||n(t)||^ + ||nt(t)||^). Here the Fourier coefficients of a 

function is given as 

= M f u{x)M^^l^dx. 


2L 


'n 


For a positive integer N the space of trigonometric polynomials of degree 
N/2 is 

Sn = span{e*''"^/^| - N/2 < k < N/2 - 1} (2.6) 

The operator —)■ Sn denotes the orthogonal projection defined as 


N/2-1 

Pnu{x) = '^ Uke 

k=-N/2 


ik'KxfL 


So that the projection operator Pn has the property 

{u — Pnu, (f) = 0, Wip E Sn- 
Pn commutes with derivative in the distributional sense: 

D^Pnu = PnD^u and D/Pnu = PnD/u. 


(2.7) 


( 2 , 8 ) 


(2,9) 
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Here and stand for the nth-order classical partial derivative with 
respect to x and t, respectively. In the rest of stndy, C denotes a generic 
constant. 


3. The semi-discrete scheme 


In this section we present a Fonrier psendo-spectral method for the spatial 
discretization and then we prove the convergence of the semi-discrete scheme. 
The equation fll.ip equivalently can be written as 

utt = ^-[u + g{u)]xx, (3.1) 

where Ti is the convolution operator given by 

nv{x)= / l3{x - y)v{y)dy. 

Jn 

Then, the equation fld.ip becomes 

= [u + g{u)]xx- (3.2) 

We dehne the operator as V. = T ^{v{k)/(]{k)) where W ^ is the 
inverse Fourier transform. The semi-discrete Fourier pseudo-spectral scheme 
for fl3.2p with the initial conditions fll.2p are given as 

n-Wt=<. + PN{9{u^)U (3.3) 

U^{x,0) = PN(j){x), uf^{x,0) = PNlpix) (3.4) 


where u^{x,t) is a function from [0,T] to Sn- We use the following lemmas 
for the convergence proof of the semi-discrete scheme: 


Lemma 3.1. 

that 


jji] For any real 0 < /i < s, there exists a constant C such 


\u - Pnu\\^ < Cm- 


\u\ 


Wu E 


(3.6) 


Lemma 3.2. Assume that f E C'^(M), u,v E n and 

k = [s] -|- 1, where s > 0. Then we have 


||/(m) - f{v)\\s < C{M)\\u - n||, (3.6) 


if ll'^^lloo < M, ||n||oo < AI, IImIIs < M and ||n||s < M, where C{M) is a 
constant dependent on M and s. 
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Lemma 3.3. Letr > 2 andv E If the kernel function f3{x) satisfies 

the condition (O, then there exists a constant C such that 


\\v\\r/2<C\\H-^/^v\\. (3.7) 

Proof. The proof directly follows from the definition of the operator 
and the condition fll.3p . □ 

The following theorem states our main result. 

Theorem 3.4. Let s > rj^, g E C'^(R) and u(x,t) he the solution of the 
periodic initial value problem fll.ip - fll.2p satisfying u(x,t) E C^([0, T]; ifp(r2)) 
for any T > 0 and u^(x,t) be the solution of the semi-discrete scheme 113. 31) - 
^ 3 . 4 ^ . There exists a constant C, independent of N, such that 

\\u - < C(T)N^/^-^\u\\x^ (3.8) 

for the initial data E Hf(Q). 

Proof. To estimate the term \\u — u^\\x ^^2 triangle inequality 

to write 

\\u - < \\u - Pnu\\x^^2 + \\Pnu - u^\\x^^2- (3-9) 

For the first term on the right-hand side, we use Lemma ITT] and the fact that 
the projection operator Pn commutes with the differentiation, to obtain the 
following estimate 

||(n - PNu)t(t)\\r /2 < CNP^-^\\ut(t)\U, (3.10) 

for s > r/2. Thus, we bound the first term in fl3.9p as 

\\u - Pnu\\x^^, < CNP^-^\\u\\x,. (3.11) 

To estimate the second term \\Pnu — u^\\x ^/2 the right-hand side of the 
inequality fl3.9p . we subtract the equation fl3.3p from fl3.ip and take the inner 
product with ip E Sx] 

{PL~^{u - u^)tt - (u- u^)xx - (g(u) - PNg(u^))xx, p) = 0. (3.12) 

Orthogonality of the projection operator Px gives 

{PL~^(u-Pxu)tu p) (3.13) 
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for all G S'tv, using this result together with fl2.8p the equation fl3.12p turns 
to 

[H~^{Pnu - u^)tt - {Pnu - u^)xx - {g{u) - Pn9{u^))xx, = 0. (3.14) 

We now set cp = {PNU — u^)t in fl3.14l) . Since is a self-adjoint operator, 

we have 

{'H~^{Pnu - u^)tt, {Pnu - u^)t) = (3-15) 

Spatial periodicity and the integration by parts yield 

{{Pnu-u^)xx, {Pnu - u^)t) = ~^\\{Pnu - u^)x{t)f- (3.16) 

For the nonhomogeneous term in fl3.14p we use the Cauchy-Schwarz inequality, 
triangle inequality and Lemma 13.21 Therefore we have 

((5'(w) - g{u^))xx, {Pnu - u^)t) 

<\{{g{u) - g{u^))x, {PNU-U^)xt)\ 

<\Mu)-g{u^))x\\ \\{PNU-u^)xtm 

<c(iie«-«'^)(()ii; + ii(PA,«-«y,(()ii?) 

< c (ii(« - PwurnWr /',+- «~)(()ii?/2 +ii(Pa« - f'')<(«)ll?) 

(3.17) 

for r > 2. Combining the results of fl3.15p - fl3.17p the equation fl3.14p becomes 

14 {\\n-^P(p^u - + ||(p„« - «"),(«)in 

< c (ii(« - p'^u)(tm„ +iKP^u - «)(()ii?,2+- «").(*)ii;) ■ 

(3.18) 

We now add the terms 

r/2 

{Pnu - {Pnu - u%) + {DUPnu - u^), D^^Pnu - u%) 
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to both sides of the inequality fl3.18p and use that 


{DliP^u - u«), DiiP^u - u%} = . 

and 

{Dl(P^u - m"), DUP^u - u%) < C {\mP„u - + \\Di(P^u - . 

Using all these results and Lemma 13.31 the inequality fl3.18p becomes 

< c (ii(« - p««)(()ii;/,+ii(p«« -«") wiiJ,2 + - t‘").(()iiL) 

< c (||(« - P„«)(«)IIL + IK-P^** - “"XOIIL + WH-^'HPnu - . 

Note that \\'H~^^‘^{PNU — u^)t{0)\\‘^ = ||(PivM —M^)(0)|p = 0. Using Gronwall 
Lemma followed by Lemma 13.11 we conclude that 

IIM-'-'XPa'U - «")X<)f + ll(P«« - t<'')(<)llL 

< f ||!i(t) - Pn!i(t)||J/ 2 e^<‘“P<iT 

Jo 

< C{T) f ||M(r)||^ dr (3.19) 

Jo 

for s > r/2. Thanks to Lemma [3.31 we know that 

ii(Pa,« - «y,(()ii?,2 + ii(p«« - 

< \\H-'/^(PmU - + |I(Pa-« - U%mr/2- 

Finally, we have the estimate 

\\Pnu - < C'(T)iV^/2-s ||^||^^_ ^3 20) 

The inequalities fl3.1ip and (13.201) give the estimation for (13.91) . which completes 
the proof of Theorem 13.41 □ 

Note that the result of the Theorem 13.41 coincides with the special cases of 
the equation (|L1D . the IBq and HBq equations, for which r = 2 and r = 4, 
respectively |l|, |l6| . 
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4. The fully-discrete scheme 


We solve the equation fll.ip by combining a Fourier pseudo-spectral method 
for the space component and a fourth-order Runge Kutta scheme (RK4) for 
time. If the spatial period [—L, L] is normalized to [0, 27r] using the trans¬ 
formation X = 7r(x + L)/L, the equation fll.ip with g{u) = becomes 

= + {n'^Y))xx- (4.1) 

The space interval [0, 27r] is divided into N equal subintervals with grid spac¬ 
ing AX = 271/N, where the integer N is even. The spatial grid points are 
then given by Xj = 27rj/X, j = 0,1,2,N. The time interval [0,T] is 
divided into M equal subintervals with grid spacing At = T/M. The ap¬ 
proximate solution to u^{Xj,t) is denoted by Uj(t). Applying the discrete 
Fourier transform to the equation fl4.ip we have 

{Uk)u = -(^Y) ^k[Uk + {Unk]. (4.2) 

The above equation can be written as an ordinary differential equation system 

{Uk)t = Vk (4.3) 

(Rfc)i = -{7Tk/Lf^k[Uk + (Unk]. (4.4) 

In order to handle the nonlinear term we use a pseudo-spectral approxima¬ 
tion. We use the fourth order Runge-Kutta method to solve the resulting 
ODE system fl4.3l) - fl4.4p in time. Finally, we hnd the approximate solution 
by using the inverse discrete Fourier transform. 


5. The Petviashvili’s iteration method 


One of the most effective numerical methods to generate solitary wave so¬ 
lution was h rst p roposed by V. I. Petviashvili for the Kadomtsev-Petviashvili 
equation in 1^. Since we do not have the exact form of the solitary wave 


solutions for the general nonlocal nonlinear wave equation, we use the Petvi¬ 
ashvili’s iteration method [l3, 13, 17 1 to construct the solitary wave solution 


numerically. 

The solitary wave solutions of fll.ll) are of the form u{x,t) = 0(a: — ct), 
where c is the propagation speed of a solitary wave. Substituting this solution 
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into and then integrating twice and using the asymptotic boundary 

conditions, we have 

cV = [/5* (0 +<^)]- (5.1) 

If we use the Fourier transform, 

1 POO poo 

u{x) = — / u{k)e^^^dk, u{k) = / u{x)e~'^^^dx, (5.2) 

2vr J_oo J-oo 

the eq. fl5.ip can be written as 

[c^ - ^{k)]<t){k) = /3{k)^P{k). (5.3) 

A simple iterative algorithm for numerical calculation of (j){k) for the above 
eq. can be proposed in the form 


k+i{k) = /^;,> ^(fc). (5.4) 

— p{k) 

where 4>n{k) is the Fourier transform of 4>n{x) which is the iteration of the 
numerical solution. Although there exists a fixed point (j)n{k) of the eq. fl5.ip . 
the algorithm fl5.4D diverges. To ensure the convergence, we add a stabilizing 
factor Mn as in [131]. The new algorithm for the nonlocal nonlinear wave eq. 
can be proposed in the form 


^+.(4) = M2 /W 

— p{k) 


where the stabilizing factor is 


~ d*:)) l‘i>n(k)?dk 

/TL Pk)<i^n{k)p{k)dk 


(5.5) 


(5.6) 


and 7 is a free parameter. Note that we can only construct the solitary wave 
solutions for the nonlocal nonlinear wave eq. under the assumption 


— /3(fc) 7^ 0. 


(5.7) 


for all k eW. 
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6. Numerical implementation 

In this section we present some numerical results to understand the effect 
of the kernel function /3(x) on the solution. Throughout the section we 
consider the quadratic nonlinearity, i.e. g{u) = v?. In the first numerical 
experiment, we consider the general kernel function satisfying fll.dp whose 
Fourier transform is 

^ 1 + + r/P sin(fc 2 ) ■ 

Here is a positive parameter. Note that (3{k) is a decreasing function of 
k. Therefore, the speed must be chosen as > 1 to satisfy the condition 
m- First, we construct a solitary wave solution corresponding to the above 
kernel for t] = 1 by using the Petviashvili’s iteration method. For this aim, 
we use the pseudo-spectral approximation for the eq. flS.Sp . The initial guess 
is 

(po{x) = ( 6 . 2 ) 

and we take the the speed as c = 1.08. In the left panel of Figure 1, we 
show the solitary wave profile of the eq. fll.ip on the interval [— 100 , 100 ] 
with N = 1024. In the right panel of Figure 1, we show the variation of 
the logarithm of the residual error with the number of iterations. Here the 
residual error is defined as 

Residual error = ||M0ji||oo (6.3) 

where 

M(j) = c^<p — [l3 * {(p + (p^)]. (6.4) 

We observe that the fastest convergence occurs when 7 = 2. 

To investigate the time evolution of the wave, the initial condition u{x, 0) 
is chosen as the profile constructed by the Petviashvili’s iteration method 
(left panel of Figure 1). On the other hand, we also need the initial condi¬ 
tion ut{x,0). The traveling wave solutions of the eq. fll.ip are given in the 
form u{x,t) = cp{x — ct). Therefore, the initial condition Ut{x,0) is chosen 
as Ut{x,0) = —cux:{x,0) and the x-derivative is evaluated by using discrete 
Fourier transform. In the left panel of Figure 2, we illustrate the time evolu¬ 
tion of the wave corresponding to these initial conditions. In the right panel 
of Figure 2, the variation of the change in the conserved quantity E (energy) 
with time is presented. 
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Figure 1: The numerical solitary wave profile for the speed c = 1.08 and the 
variation of the logarithm of the residual error with the number of iterations. 



Figure 2: The time evolution of the wave generated by Petviashvili’s iteration 
method (left panel) and the variation of the change in the conserved quantity E 
(energy) with time (right panel). 


We point out that the kernel (16.ip corresponds to the kernel giving rise to 
the IBq equation when rj = 0. In this case, if the nonlinearity is quadratic, 
the relation between the amplitude A and the speed parameter c is given 
hj A = 1.5(c^ — 1). In Figure 3, the dashed line shows the variation of the 
amplitude with the speed parameter for the IBq equation. The solid lines 
show the variation of the amplitude with the speed parameter corresponding 
to kernel (16.Ih for rj = 10, rj = 5, rj = 1 and rj = 0.1. 
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Figure 3: Variation of the amplitude with the speed parameter for the IBq eq. 
and for the nonlocal nonlinear wave eq. with the kernel (I6.ip for various values of 
rj. 


The question arises naturally how the solitary wave solutions of the nonlocal 
nonlinear wave equation behaves when 77 —)■ O’*". For this aim, we perform 
some numerical tests for various values of rj by using the initial condition 


u{x, 0) 
v{x, 0) 


-sech^f 

4 

- sech^ 
2 






(6.5) 

( 6 . 6 ) 


The given initial condition corresponds to the exact solitary wave solution of 
the IBq equation. The problem is solved on the interval —100 < x < 100 for 
times up to T = 10. In Figure 4, the solid line shows the exact solution of the 
IBq equation given in [l| which corresponds to the initial conditions fl6.5p - 
fl6.6p . The dashed lines show the numerical solution of the nonlocal nonlinear 
wave equation obtained by the kernel function whose Fourier transform is 
given by fl6.1l) for t] = 10, rj = 5, rj = 1 and rj = 0.1 at the final time 
T = 10. The numerical tests show that the solutions of the nonlocal nonlinear 
wave equation converge to the solitary wave solution of the IBq equation as 
r] —)■ O’*". 
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Figure 4: The exact solution of the IBq equation and the numerical solution of the 
nonlocal nonlinear wave equation with ()6.ip . 


Next, we consider the kernel function whose Fourier transform is 




1 yU 

l + k? + k‘ ^ 1 + k*' 


(6.7) 


Here /r is a positive parameter. As (3{k) is a decreasing function of k, the 
speed must be chosen as > 1 + /i to satisfy the condition fl5.7p . This kernel 
corresponds to the kernel giving rise to the HBq equation with rji = ri 2 = 1 
when /i = 0. To observe the behavior of the solutions as /i —?• O’*", we perform 
some numerical tests for various values of fi by using the initial condition 

1 T 

u{x,0) = -sech^f—=), (6.8) 

266 ^2x/l3 

<'(-.0) = ( 5 ^ ) tauh ) . (6.9) 

The given initial condition gives rise to the exact solitary wave solution of the 
HBq equation. In Figure 5, the solid line shows the exact solution of the HBq 
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equation given in which corresponds to the initial conditions fl6.8p - fl6.9p . 
The dashed lines show the numerical solution of the nonlocal nonlinear wave 
equation obtained by the kernel function whose Fourier transform is given 
in fl6.7p for p = 10, /i = 5, /i = 1 and /i = 0.1 at the final time T = 10. 
The numerical tests show that the solutions of the nonlocal nonlinear wave 
equation converge to the the solution of the HBq equation as /r —)■ O'*". 




Figure 5: The exact solution of the HBq equation and the numerical solution of 
the nonlocal nonlinear wave equation with (16.71) . 

For the last experiment, our aim is to investigate whether the above 
convergence property is also valid for the blow-up solutions or not. The 
blow-up theorem for the nonlocal nonlinear wave equation given in has 
been restated as follows: 

Theorem 6.1. (Theorems 3.4 and 5.2 off^]) Let s > 1/2 and r >2. Then, 
there is some T > 0 such that the Cauchy problem fll.l|) - fll.2|) is well-posed 
with solution in ^^([0, T], for initial data </>, t/ G 77^(12). Further¬ 
more, suppose that P(j),P'ilj G and G 77^(12). If there is some 
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1 / > 0 such that 


pg (yP) < 2z/p^ + 2 (1 + 2z/) G {p) for all p G M, 
and the initial energy 

E (0) = \\Pfjf + Uf + 2 [ Gi(l))dx<0 , 

Jm. 

then the solution u of the Cauchy problem blows up in finite time. 

For the numerical tests, we use the kernel whose Fourier transform is given 
in flh.ip and we choose the initial conditions 

/ \ 2 2 
(f){x) = 4 f — - 1 j e 3 “, ^p{x) = {x^ — l) e 3“ ( 6 . 10 ) 

as in j^. Thus, the conditions for the blow-up theorem are satished for 
z/ = 4. The problem is solved on the interval —10 < a; < 10 for times up 
to T = 3. In Figure 6 , we present the variation of the Loo norm of the 
numerical solution obtained by using the Fourier pseudo-spectral scheme. 
The numerical results strongly indicate that the blow-up time of the nonlocal 
nonlinear wave equation converges to t* ~ 1.8, which is the blow-up time of 
the IBq equation, as p —)■ O'*". 
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Figure 6: Loo norm of the numerical solution for the increasing time. 
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